Landau level degeneracy and quantum Hall effect in a graphite bilayer 
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We derive an effective two-dimensional Hamiltonian to describe the low energy electronic exci- 
tations of a graphite bilayer, which correspond to chiral quasiparticles with a parabolic dispersion 
exhibiting Berry phase 2tt. Its high-magnetic-field Landau level spectrum consists of almost equidis- 
tant groups of four-fold degenerate states at finite energy and eight zero-energy states. This can be 
translated into the Hall conductivity dependence on carrier density, a xy (N), which exhibits plateaus 
at integer values of 4e 2 //i and has a "double" 8e 2 /h step between the hole and electron gases across 
zero density, in contrast to {An + 2)e 2 /h sequencing in a monolayer. 



PACS numbers: 73.63.Bd, 71.70.Di, 73.43.Cd, 81.05. Uv 
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For many decades, the electronic properties of a 
graphite monolayer have attracted theoretical interest 
due to a Dirac-type spectrum of charge carriers 0, 0, 
IE QL Ej hi this gapless semiconductor 0. Recently 
Novoselov et al. [8j fabricated ultra-thin graphitic de- 
vices including monolayer structures. This was followed 
by further observations 0, IToL ITU of the classical and 
quantum Hall effects (QHE) in such systems confirming 
the expectations Q of an unusual phase of Shubnikov 
de Haas oscillations and QHE plateaus sequencing, as 
manifestations of a peculiar magneto-spectrum of chiral 
Dirac-type quasiparticles containing a Landau level at 
zero energy 

In this Letter we show that quasiparticles in a graphite 
bilayer display even more intriguing properties includ- 
ing a peculiar Landau level (LL) spectrum: these are 
chiral quasiparticles exhibiting Berry phase 2n, with a 
dominantly parabolic dispersion and a double-degenerate 
zero-energy LL incorporating two different orbital states 
with the same energy. Taking into account spin and val- 
ley degeneracies, the zero-energy LL in a bilayer is 8- 
fold degenerate, as compared to the 4-fold degeneracy 
of other bilayer states and the 4-fold degeneracy of all 
LLs in a monolayer. The structure and degeneracies of 
the Landau level spectrum in a bilayer determine a spe- 
cific sequencing of plateaus in the density dependence 
of the QHE conductivity a xy (N) which is distinguish- 
ably different from that of Dirac-type quasiparticles in a 
graphite monolayer and of non-chiral carriers in conven- 
tional semiconductor structures. 

We model a graphite bilayer as two coupled hexagonal 
lattices including inequivalent sites A, B and A, B in the 
bottom and top layers, respectively. These are arranged 
according to Bernal (A-B) stacking 0, fl^L Il3| , as shown 
in Fig.n A lattice with such symmetry supports a degen- 
eracy point at each of two inequivalent corners, K and K, 
of the hexagonal Brillouin zone 01 > which coincide with 
the Fermi point in a neutral structure and determine the 
centers of two valleys of a gapless spectrum. At the de- 
generacy point, electron states on inequivalent (A/ B or 
A/B) sublattices in a single layer are decoupled, whereas 
interlayer coupling j^b = T 1 f° rms 'dimers' from pairs 
of A-B orbitals in a bilayer [solid circles in Fig. thus 




FIG. 1: Left: schematic of the bilayer lattice (bonds in the 
bottom layer A, B are indicated by solid lines and in the top 
layer A, B by dashed lines) containing four sites in the unit 
cell: A (white circles), B (hashed), AB dimer (solid). Right: 
the lattice of a monolayer. 



leading to the formation of high energy bands [T3. fl3| . 
The low energy states of electrons are described by 
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The effective Hamiltonian H2 operates in the space of 
two-component wave functions <!> describing electronic 
amplitudes on A and B sites and it is applicable within 
the energy range |e| < j7i- In the valley K, £ = +1, 
we determine $^=+i = (<j)(A), 4>(B)), whereas in the val- 
ley K, £ = —1 and the order of components is reversed, 
<f>£ = _i = (4>(B),(fi(A)). Here, we take into account two 
possible ways of A ^ B hopping: via the dimer state 
(the main part) or due to a weak direct AB coupling, 
7 Ai j = 73 <C (the term h w ). They determine the 
mass m — 7i/2v 2 and velocity to = (v3/2) a/y A jj/h. 
Other weaker tunneling processes pj are neglected. The 
term h a takes into account a possible asymmetry between 
top and bottom layers (thus opening a mini-gap ~ u). 
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For comparison, the monolayer Hamiltonian 0], 

Hi = £V ^ ^ EE £ V + (JyPy) , 

is dominated by nearest neighbor intralayer hopping 
Iab = 1b A ee 7q > 7x, so that v = (V3/2) ajAB/h. 
For equivalent parameters in bulk graphite 13, «3 <C v 
in Eq. Thus, the linear term h w , which is similar 
to Hi, is relevant only for very small electron momenta 
(i.e., in an electron gas with a small density at a very 
low magnetic field) whereas the energy spectrum within 
the interval 2~Ji(v3/v) 2 < \e\ < jji is dominated by the 
first term |15| in H 2 producing a dispersion |e| = p 2 /2m, 
which contrasts with |e| = vp in a monolayer. Hi and 

i?2 form a family of Hamiltonians Hj — £, J f(\p\) cr-n de- 
scribing particles which are chiral in the sublattice space, 
where n = 1^ cos(Jtp) +l y sin(J^) for p/p = (cos <p, sin <p>) 
[tt = pe lcp ], and the degree of chirality is J = 1 in a 
monolayer and J = 2 in the bilayer. It is interesting to 
notice that quasiparticles described by the Hamiltonians 
Hj acquire a Berry phase Jtt upon an adiabatic prop- 
agation along a closed orbit, thus charge carriers in a 
bilayer are Berry phase 2tt quasiparticles, in contrast to 
Berry phase tt particles in the monolayer of graphene 0] . 
According to the inverted definition of sublattice compo- 
nents (for which we reserve 2x2 Pauli matrices oi [lfil ]) 
of the wave functions $j = +i and quasiparticles 
in different valleys, £ = ±1 have effectively the oppo- 
site chirality. Also, the existence of two valleys is crucial 
for the time-reversal symmetry of the chiral Hamiltoni- 
ans. In application to H2, time reversal is described by 
(Hi <g) a x )H* (p, B, u) {Hi ®a x ) = H (-p, -B^u), where 
Hi swaps £ = +1 and £ = — 1 in valley space |16| . 

A microscopic analysis leading to the bilayer Hamilto- 
nian H2 uses the tight-binding model of graphite and the 
Slonczewski-Weiss-McClure parameterization |7j] of rele- 
vant couplings. We represent the Hamiltonian near the 
centers of the valleys in a basis with components corre- 
sponding to atomic sites A, B, A, B in the valley K [blj 
and to B,A,B,A in the valley K, and distinguish be- 
tween on-site energies, ±5" in the two layers, 
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\ [ K=p x +lPy, 

p = — ihV — eA, 
B = rotA, 
J { = 2heB. 



The Hamiltonian TL determines the following spectrum 
of electrons in a bilayer at zero magnetic field. There are 
four valley-degenerate bands, e^(p), a = 1,2, with 
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The dispersion ei(p) describes low energy bands. In 
the intermediate energy range, ^71 (^3/^) ,u < \ei\ < 
71 , it can be approximated with 
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This corresponds to the effective mass for electrons near 
the Fermi energy in a 2 D gas with dens ity N, m c — 
p/{dei/dp) = (7i/2w 2 ) y/l + 4nh 2 v 2 N/j 2 . The relation 
in Eq. (|5J) interpolates between a linear spectrum t\ « vp 
at high momenta and a quadratic spectrum t\ ss p 2 /2m, 
where m = 71 /2u 2 . Such a crossover happens at p ~ 
ji/2v, which corresponds to the carrier density ./V* ss 
7i/ (47r/i 2 7j 2 ). The experimental graphite values 0. ITo| 
give N* w 4.36 x 10 12 cm~ 2 , whereas the dimer band 
62 becomes occupied only if the carrier density exceeds 
Ar( 2 ) « 2 7l 2 / (irh 2 v 2 ) w 8N* w 3.49 x 10 13 cm- 2 . The 
estimated effective mass m is light: m m 0.054m e using 
the bulk graphite values El ■ 

The 4x4 Hamiltonian TL contains information about 
the higher energy band £2, and, therefore, is not conve- 
nient for the analysis of transport properties of a bilayer 
which arc formed by carriers in the low energy band ei. 
We separate TL into 2x2 blocks, where the upper left di- 
agonal block is Hn ee £,(^ua z +v 3 [(T x p x —a y Py}), the lower 
right diagonal block is H22 = — ^£,u(T z +'-fia x , and the off- 
diagonal blocks are H21 = H12 = v^(a x p x + a y p y ). Then, 
we take the 4x4 Green function determined by TL, eval- 
uate the block Gn related to the lower-band states, and 
use it to identify the effective low-energy bilayer Hamil- 
tonian i?2- Using Gal = {H aa — e) 1 , we write 
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Then, we find that Gn = (l - G^ Hi 2 G^ H 2 i\ G$ , 

so that Gil + £ = Hn — H^G^ H 2 i- Since |e| -C 71, we 

expand G 2 °2 = (H22 — e) 1 in 7] -1 , keeping only terms 
up to quadratic in p (and therefore in ir^ir), and arrive 
at the expression in Eq. QJ. 

For low quasiparticlc energies, |e| <C 71, the spec- 
trum determined by H 2 in Eq. agrees with ei (p) 
found using the 4x4 Hamiltonian TL. Similarly to bulk 
graphite 0,0, the effect of h w consists of trigonal warp- 
ing, which deforms the isoenergetic lines along the direc- 
tions ip — ipo, as shown in Fig. [21 For the valley K, 
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where €2 describes the higher-energy (AB dimer) bands. 



|7r. At the lowest energies |e| < ^71 (v 3 /v) , trigonal 
warping breaks the isoenergetic line into four pockets, 
which can be referred to as one "central" and three "leg" 
parts 0| . The central part and leg parts have minimum 
|e| = \u at p — and at \p\ — 71W3/W 2 , angle (fo, re- 
spectively. For i> 3 ~ 0.1c. wc find (using the data in 
Ref. 0) that the separation of a 2D Fermi line into four 
pockets would take place for very small carrier densities 
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FIG. 2: Schematic of the Fermi line in the valley K , £ = 1, 
for high (light shading) and low density (dark shading). Note 
that the asymmetry of the Fermi line at valley K, £ = —1, is 
inverted. 



N < N c = 2(v 3 /v) 2 N* ~ 1 x 10 n cTO- 2 . For N < N c , 
the central part of the Fermi surface is approximately 
circular with area A c ~ tvs 2 /(hv 3 ) 2 , and each leg part is 
elliptical with area A\ « \A C - This determines the fol- 
lowing sequencing of the first few LL's in a low magnetic 
field, B < B c w hN c /Ae ~ IT. Every third Landau level 
from the central part has the same energy as levels from 
each of the leg pockets, resulting in groups of four degen- 
erate states. These groups of four would be separated by 
two non-degenerate LLs arising from the central pocket. 

In structures with densities N > N c or for strong mag- 
netic fields B > B c , the above described LL spectrum 
evolves into an almost equidistant staircase of levels. We 
derive such a spectrum numerically from Eq. Q using the 
Landau gauge A = (0, Bx), in which operators 7r and it 
coincide with raising and lowering operators [Tsl ] in the 
basis of Landau functions e lky (j) n (x), such that n^(f> n — 
i(h/X B ) y/2(n + l)^„+i, 7T0„ = -i(h/X B )V2n4> n -i, and 
tv (fio =0, where X B — yh/JeB). In this we followed an 
approach applied earlier to bulk graphite yjfl and used 
the bulk parameters for intralayer v and interlayer 71, 
varying the value of the least known parameter v 3 . The 
spectrum for the valley K (£ = 1) is shown in Fig. |3 
as a function of the ratio v 3 /v for two different fields. 
Fig. OJa) shows the evolution of the twenty lowest levels 
for B = 0.1T as a function of v 3 , illustrating the above- 
mentioned crossover from an equidistant ladder at v 3 = 
to groups of pocket-related levels. 

The LL spectrum obtained for B = IT, Fig. G^b) re- 
mains independent of v 3 over a broad range of its val- 
ues. Hence, even in the absence of a definite value of 
V3, we are confident that the LL spectrum in bilayers 
studied over the field range where ^A^ 1 > v 3 m can be 
adequately described by neglecting v 3 , thus using an ap- 
proximate Hamiltonian given by the first term in H 2 , 
Eq. Q ■ The resulting spectrum contains almost equidis- 
tant energy levels which are weakly split in valleys K 
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FIG. 3: Numerically calculated Landau levels for various val- 
ues of U3 for valley K, £ = 1. (a) for B = 0.1T (huj c = 
0.216meV assuming 71 = 0.39el/, v = 8.0 x 10 5 m/s), (b) 
for B — IT (hu} c — 2.16meV). Broken lines show groups of 
four consecutive levels that become degenerate at large V3, 
with each group separated from the next by two solid lines 
representing levels that are not degenerate. 

(£ = +1) and K (£ = -1), 

et = ±hu c ^n(n- 1) - ±£<5, for n > 2, (3) 

Here, w c = eB/m, e+ and e~ are assigned to electron 
and hole states, respectively, and D n £ = [e — £,u/2 + 
£n5]/(hu c y/n{n- 1)), C ni = l/^l + \D ni \ 2 . In the 
limit of valley (u — 0) and spin degeneracies [13, we 
shall refer to these states as 4-fold degenerate LLs. 

The LL spectrum in each valley also contains two levels 
identified using the fact that 7r 2 0i = Tv 2 <fio = 0, 

Je = |^u; $ O £ = (0o,O); 

\ ei = |£«-#; ^ = (^,0). w 

According to different definitions of two-component $ 
in two valleys, n = 0, 1 LL states in the valley K are 
formed by orbitals predominantly on the A sites from 
the bottom layer, whereas the corresponding states in 
the valley K are located on B sites from the top layer, 
which is reflected by the splitting u between the lowest 
LL in the two valleys. In a symmetric bilayer (u = 0) 
levels £0 and e\ are degenerate and have the same energy 
in valleys K and K, thus forming an 8-fold degenerate 
LL at e — (here, spin is taken into account). Also, 
note that the spectrum of high-energy LLs, Eq. © is 
applicable in such fields that hX^ 1 < ji/2v. For higher 
fields the full two-band Hamiltonian Ti has to be used to 
determine the exact LL spectrum, nevertheless, the 8-fold 
degeneracy of the zero-energy LL remains unchanged. 

The group of 8 states at |e| = (4 for electrons and 
4 for holes, Eq.QJ) embedded into the ladder of 4- fold 
degenerate LL's with n > 2, Eq. © is specific to the 
magneto-spectrum of J = 2 chiral quasiparticles. It 
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FIG. 4: Landau levels for a bilayer (left) and monolayer 
(right). Brackets (n,£) indicate LL number n and valley in- 
dex £ = ±1. In the center the predicted Hall conductivity 
(7 X y (center) as a function of carrier density for bilayer (solid 
line) is compared to that of a monolayer (dashed line). 



would be reflected by the Hall conductivity dependence 
on carrier density, <7 xy (N) shown in Fig. A solid 

line sketches the form of the QHE a xy ' (N) in a bilayer 
which exhibits plateaus at integer values of 4e 2 /h and 
has a "double" 8e 2 /h step between the hole and elec- 
tron gases across N = that would be accompanied 
by a maximum in a xx ■ Figure 0] is sketched assuming 
that temperature and the LL broadening hinder small 
valley and spin splittings as well as the splitting be- 
tween n = 0,1 electron/hole LL's in Eqs. QJ, so that 
the percolating states from these levels would not 
be resolved. To compare, a monolayer has a spectrum 



containing 4- fold (spin and valley) degenerate LLs pj, 
eo = and £^>i — ±V%nfiv / \b shown on the r.h.s of 
Fig. ^] which corresponds to Hall conductivity a x ]J (TV) 
exhibiting plateaus at (An + 2)e 2 Jh (dotted line |2(J), as 
discussed in earlier publications [3j. 

The absence of a a xy = plateau in the QHE accom- 
panied by the maximum in a xx in the vicinity of zero 
density is the result of the existence of the zero-energy 
LL, which is the fingerprint of a chiral nature of two- 
dimensional quasiparticles. This contrasts with a grad- 
ual freeze-out of both Hall and dissipative conductivities 
in semicondutor structures upon their depletion. Hav- 
ing compared various types of density dependent Hall 
conductivity, we suggest that two kinds of chiral (Berry 
phase Jtt) quasiparticles specific to monolayer (J = 1) 
and bilayer (J = 2) systems can be distinguished on the 
basis of QHE measurements. It is interesting to note 
that the recent Hall effect study of ultra-thin films by 
Novoselov et al. [l(| featured both types of a xy (N) de- 
pendence shown in Fig. 

It is also worth mentioning that the 8-fold degeneracy 
of the group of e = LLs in a bilayer, Eqs. I@J), is quite 
unusual in 2D systems. It suggests that e-e interaction in 
a bilayer may give rise to a variety of strongly correlated 
QHE states. For structures studied in Ref. flfl with elec- 
tron/hole densities N ~ 10 12 cm~ 2 , such a regime may 
be realized in fields B ~ 10T. 
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